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1 Introduction 


Precise measurements of the parameters of the Cabibbo-Kobayashi-Maskawa unitarity 
triangle are of great value in searching for manifestations of new physics in the flavour 
sector. In particular, the determination of the angle 7 = arg(—14dKib/^dKb) (^^^o denoted 
as 03 in the literature) in processes involving tree-level decays provides a Standard Model 
(SM) benchmark against which observables more sensitive to new physics contributions can 
be compared. Currently such comparisons are limited by the uncertainty on 7 , which is 
~ 7° [2]-[^. More precise measurements and new analysis strategies are therefore required. 

Sensitivity to 7 in tree-level processes may be obtained through the study of CP- 
violating observables in the decays Dh^, where D indicates a neutral charm meson 

which decays in a mode common to both and states, and h, the bachelor hadron, 
is either a kaon or a pion. In the case of B~ — DK~, interference occurs between the 
suppressed b —)■ ucs and favoured b —)■ cus quark-level transitions, and similarly for the 
charge-conjugate decay. The magnitude of the interference is governed by three parameters: 
the weak-phase difference, 7 , the CP-conserving strong-phase difference, 6b, and the ratio 
of the magnitudes of the two amplitudes, r^. Similar interference effects occur in the 
case when the bachelor hadron is a pion, but then additional Cabibbo suppression factors 
mean that the sensitivity to 7 is much reduced. Many classes of D decay can be exploited. 
Important examples include the so-called ADS modes j^, which are decays to quasi 
flavour eigenstates such as D —)■ and the GLW modes (^, which are decays to CP 

eigenstates such as D —?■ K^K~. Measurements exist from LHCb that follow both the 
ADS and GLW approaches [8 
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as well as alternative methods 12,13 


In the case that the D meson decays to three or more hadrons, the interference effects 
that are sensitive to 7 vary over the phase space of the D decay due to the role of 
strongly-decaying intermediate resonances. If the D decay is analysed inclusively, the 
integration over phase space in general dilutes the net sensitivity. For multibody ADS 
modes the dilution factor can be measured with DD pairs coherently produced at the 

LHGb has previously made use of such measurements performed 

in B^ — Dh^ analyses exploiting the 


0(3770) resonance 14 


with data from the GLEO-c experiment 15 


TT 


17 

T 


10 . It has recently been pointed out 18 


modes D —and D — K'^K^ 
that similar considerations apply to self-conjugate multibody modes such as D ^ 

These modes approximate to CP eigenstates, and hence a B^ — DK^ analysis that 
employs them can be considered a quasi-GLW (qGLW) analysis. In this case the dilution 
factor is related to how closely the mode approaches a CP eigenstate, and can also be 
measured at the open charm threshold. 

This paper presents the measurement of CP observables from B^ —)■ Dh^ decays, 
where D mesons are reconstructed using three different multibody hnal states. These 
decays are the ADS channel D —)■ and the quasi-GLW modes D —)• 7 r''' 7 r“ 7 r‘^ and 

D A+A-ttF In all cases, higher sensitivity is attained compared with the results of 


the previous measurements which exist from the BaBar 19 and Belle 20 collaborations 
for the ADS channel, and from BaBar for the mode D —)■ tt’ 


- — 0 
TT TT^ 


21 . Measurements at 


the 0(3770) resonance indicate that the dilution effects in D —)■ and 
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D —)■ 7T~^'n'~7T^ are rather small, making these decays particularly suitable for an inclusive 
analysis. 

This paper is organised as follows. Section introduces the observables that the 
analysis seeks to measure, and explains how they are related to the underlying physics 
parameters and the dilution factors that are determined externally. Section describes the 
LHCb detector and the data set on which the analysis is based. Sections and present 
the candidate selection and the analysis procedure. Results are given in Sect. together 
with a discussion of the systematic uncertainties. In Sect, [^the measured observables are 
interpreted in terms of 7 and the other physics parameters, and conclusions are drawn. 


2 Observables and external inputs 


In the ADS channel there exist two suppressed modes, —>■ and two 

favoured modes, —)■ [K^n^'n'^]Dh^ , for h = K and vr. In both cases the suppressed 

modes are as yet unobserved, although Belle has reported hrst evidence for B^ —>■ 
[71^K^7i^]j:iK^ and —>■ [Tr^A'^Tr^jovr^ |^. As is customary in an ADS analysis, the 

“““ _ r(i;- ^ ln-K+7r'‘]Dh-) + r(i;+ ^ ln+K-n«]oh+) 

ADS(l) [x-n+7r0]oh-) + r(B+ -» [K+7r-n<‘]oh+) ^ ’ 

is dehned to give the relative rates of the suppressed to the favoured decays. The asymmetry 

,K„o ^ r(B- ^ - r(g+ ^ lT+K-T«]oh+) 

ADS(k) ^ [r-K+7r<>]nh-) + r{B+ -» [ir+A'-ir»]c/i+) ^ 

quantihes the amount of CP violation in the suppressed modes. An asymmetry is also 
constructed for the favoured channels. 


^ T{B- ^ [K-TT+Tr^jpK-) - r(R+ ^ [i^+7r-7r>ir+) 
^ “ T{B- [K-7r+7rO]j:>K-) + T{B+ [K+7r-7rO]DK+)' 


The observables ^ads(x) carry the highest sensitivity to the angle 7 ; they 

depend on the underlying physics parameters as 


^Ar>S{K) 

A i^TTTT*^ 

^ADS(i4:) 






cos(<5s + ^d""°)cos7. 


sin(5B + sin 7 


/ TyKirTT^ 

/ -^ADS(A')* 


( 4 ) 

( 5 ) 


Here ~ 0.05 is the ratio of the magnitudes of the doubly Cabibbo-suppressed and 
Cabibbo-favoured D decay amplitudes and 6^'^'^° is the strong-phase difference between 
the amplitudes, averaged over phase space. The coherence factor accounts for 

possible dilution effects in the interference arising from the contribution of the intermediate 
resonances in the D decay . Both and have been measured with quantum- 

correlated DD decays collected at the '0(3770) resonance by the CLEO-c experiment. 
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, where the phase- 


and have been found to be and 0.82 ± 0.07, respectively 

difference is given in the convention where CP\D^) = \D^). The relatively large 

value of means that the dilution effects are small, and hence this decay is a promising 

mode to exploit for the measurement of 7 . Note that for reasons of clarity Eqs. and 
are restricted to terms of O ({tbY, small effects of 


mixing are omitted. Full expressions may be found in Ref. 22 


In the quasi-GLW analysis of the two self-conjugate modes D — h'^h'~TT^ [h' = K,7r), 
observables are dehned analogously to those used in the CP-eigenstate case. The hrst of 
these is the ratio of partial widths 


R 


T{B- 




■) + r(p+ ^ D 

h' K+) 


qGLW 


r{B- D^K-) + r{B+ D^K+) 


( 6 ) 


where D^h'h'nO signihes a D meson with fractional CP-even content P: 


h'h'7T° 


Both the 

numerator and the denominator of Eq. involve B meson partial widths only and have 
no dependence on the D meson branching fractions. In practice, therefore, PqGLw 
determined by forming the ratio of two more ratios. 


P: 


R 


h'h'w° 

qGLW 

h'h'iT^ 

Kj-K 


R 


Kj-n 


^ / R 

! R 


Kj-n 


r(p- ^ [pp7r%p-) + r(p+ ^ [hpv%p+) 


R 


Kirn^ _ 

K jn 


r(P- ^ [PPvrOjzjTT-) + r(P+ ^ [PP7r0]fl7r+) ’ 
r(p- ^ [p-7r+7r%p-) + r(P+ ^ [K+n-Tr^]DK+) 
T{B~ —)■ [K~Tr+Tr^]DTr~) + T{B+ —)■ [P+7r“7r0]£)7r+) 


(7) 

( 8 ) 

(9) 


where the approximate equality in Eq. acknowledges that very small interference effects 
in the B^ —)■ decays specihed in Eqs. and can be neglected. This is a good 

assumption because the ratio between interfering amplitudes in B^ — Dtt^ decays is 
known to be very small 4 . Furthermore, the ratio R^c^Y^ / 7*6 interpreted in 
terms of the underlying jmysics parameters, taking these interference effects into account. 
Asymmetries, PqGLw(h) ~ ^ire also constructed, where 


A 


h'h'nO _ r(p- ^ [Cp7r%h-) - r(P+ ^ [pp7r%h+) 


qGLW(/i) 


( 10 ) 


r(P- ^ [PCvrOj^h-) + r(P+ ^ lh'h'7r0]Bh+)' 

The relations between RYYini and the most sensitive to 7 of the two 


. ^ and Ah'h'w^ 

''qGLW ^qGLW(i<') 

asymmetries, and the underlying physics parameters are 


-Qh'h'TT^ 

■'^qGLW 

A h'h''K^ 

^qGLW(iC) 


= 1 + {tbY + (2P 

?h'h'-K° 


h'h'-IT° 


= (27 


l)2rB cos Sb COS 7 , 
l) 2 rB sin 5b sin 7 /PJglw- 


( 11 ) 

( 12 ) 


The small effects of mixing are neglected, but can be accommodated if required 

A recent analysis using CLEO-c data 


18 


18 


has used decays of coherently produced DD 


pairs to determine Pf 


= 0.968 ± 0.018 and F 


K+K- 


= 0.731 ± 0.062. The high 
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value of F' 


'TT^TT TT 


implies that the decay —)■ tt 


- — 0 • 
TT TT^ ^ 


is very close to being a CP-even 


eigenstate and the interference terms in Eqs. 11 and 12 suffer very little dilution, tending 


towards the equivalent GLW CP-even expressions. 

When measuring CP asymmetries at the LHC, it is necessary to allow for the possibility 
that the initial state may contain different numbers of B~ and mesons. Therefore, a 
production asymmetry, 

- a{B-) + a(B*Y ' 

is dehned where <j{B~) and <j{B^) are the cross-sections for the production of B~ and B^ 
mesons, respectively, within the LHCb acceptance. 

To summarise, twelve observables are measured in total: the two ADS asymmetries 
^ADS(h)’ ratios Plnsw asymmetry the four quasi-GLW asymmetries 

^qGLW(/i) ratios PqCLw! B^/B~ production asymmetry, Aprod- 


3 The LHCb detector and data set 

The analysis uses data collected by LHCb in pp collisions at ^/s = 7 TeV in 2011 and 8 TeV 
in 2012, corresponding to integrated luminosities of 1.0 fb“^ and 2.0 fb“^, respectively. The 
LHCb detector 23,24 is a single-arm forward spectrometer covering the pseudorapidity 


range 2 < p < 5, designed for the study of particles containing 6 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The polarity of the dipole 
magnet is reversed periodically throughout data-taking in order to combat systematic 
biases due to possible detector asymmetries. The tracking system provides a measurement 
of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% 
at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary 
vertex, the impact parameter, is measured with a resolution of (15 -|- 29/pt) M-ni, where px 
is the component of the momentum transverse to the beam, in GeV/c. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors. Photons, electrons and hadrons are identihed by a calorimeter system consisting 
of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identihed by a system composed of alternating layers of iron and 
multiwire proportional chambers. The online event selection is performed by a trigger 25 


which consists of a hardware stage, based on information from the calorimeter and muon 
systems, followed by a software stage, which applies a full event reconstruction. Offline 
a loose selection based on a decision tree algorithm 26 is run to reduce the size of the 
sample prior to final analysis. 

Approximately one million simulated events (after geometric detector acceptance) of 
each class of signal decay are used in the analysis, as well as a large inclusive sample 
of generic Bg —)■ DX decays, where q G {u,d,s}. In the simulation, pp collisions are 
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generated using P ythia 27 with a specific LHCb configuration . Decays of hadronic 

in which hnal-state radiation is generated using 


particles are described by EvtGen 
Photos 
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30 . The interaction of the generated particles with the detector, and its response, 


are implemented using the Geant4 toolkit 31 as described in Ref. 32 


4 Candidate selection 


The events used in the analysis must be selected by the hardware trigger, either for the 
case where the candidate triggered the event via the hadronic calorimeter (and not 
the muon system), or where the event was triggered independently of the candidate. 
The study is performed with B^ — Dh^ candidates, where the D meson decays to a 
three-body hnal state composed of any combination of two charged kaons and pions 
and a 7r° candidate. The 7r° is identihed by a decay to two photons, as recorded by the 
electromagnetic calorimeter. 

All candidates passing the B^ Dh^ reconstruction are required to have an invariant 
mass in the range of 5080-5900 MeV/c^. The mass of the reconstructed D candidate is 
required to be within ±50 MeV/c^ of the nominal mass 1^. In addition, the mass of 


the 7r° candidate must be within ±20 MeV/c^ of the nominal vr^ mass j^. Both of these 
mass windows correspond to approximately plus or minus twice the mass resolution of 
the respective reconstructed particles. The 7r° candidate must also have a momentum 
of pt > 0.5 GeV/c and p > 1.0 GeV/c. The bachelor particle is required to satisfy 
0.5 < pt < 10 GeV/c and 5 < p < 100 GeV/c, while the charged D daughters must have 
Pt > 0.25 GeV/c. In order to improve the resolution of the mass of the B^ candidate, the 


decay chain is rehtted 33 constraining the positions of the B^ and D vertices, while at 


the same time constraining the D candidate to its nominal mass. 

In addition to these selection criteria, further background suppression is achieved 


through the use of a boosted decision tree (BDT) discriminator 34 using the GradientBoost 
algorithm 35 . The BDT is trained using a signal sample of B^ —events from 


simulation and a sample of pure combinatorial background from data with B^ candidates’ 
invariant mass greater than 5900 MeV/c^, which are not used in the invariant mass ht. The 
BDT utilises a variety of properties associated to each signal candidate. These properties 
include: p and px of the D meson, the D daughter candidates and the bachelor particle; 
and the Xip of fhe D meson, charged D daughter candidates, bachelor particle and the B^ 
meson (where Xip is dehned as the difference between the of the primary vertex (PV) 
reconstructed with and without the particle of interest). Other properties include: the 
flight distance from the PV for the B^ and D candidates; vertex quality, por degree 
of freedom, for the B^ and D candidates; and the angle between the line connecting the 
PV to the particle’s decay vertex and the particle’s momentum vector for the B^ and D 
candidates. Another characteristic used in the BDT is an isolation variable representative 
of the Pt imbalance surrounding a B^ candidate. The variable is dehned as 




Pt(B^) - EnPT 




(14) 


5 











where the sum is performed over the n tracks lying within a cone around the candidate, 
excluding the tracks related to the signal. The cone is dehned by a circle of radius f.5 
units in the plane of pseudorapidity and azimuthal angle (measured in radians). No PID 
information is used as an input variable; consequently the BDT has similar performance 
for both the —)■ DK^ and —)■ decay modes, with some slight variation arising 

due to differences in kinematics between the two. 

The optimal cut value of the BDT is determined by optimising the metric s/\/s + b, 
where s is the expected signal yield in the suppressed B^ —)■ DK^ ADS mode and b is 
the combinatoric background level as taken from the favoured mode, which is expected to 
have comparable background levels to the snppressed mode. The expected signal yield 
is calcnlated as the yield in the favoured B^ —)■ Dn^ ADS mode scaled by the predicted 
branching fraction of the B^ —)■ DK^ mode and by the expected ratio between the 
suppressed and favoured ADS modes, while taking into account differences in PID efficiency. 
Assessment of this s/y/s + b metric finds a working point where a signal efficiency of 85% 
is expected while rejecting > 99% of combinatorial backgronnd. A similar optimisation 
procednre performed nsing the 5^—)- [7r~^7r~7r^]DK^ and 5^—)- [K~^K~n^]£)K^ decays 
retnrns a comparable working point, and thns the same reqnirement is imposed in the 
selection of the quasi-GLW modes, as well as the ADS modes. 

Particle identification, essential for the distinction between B^ —> DK^ and B^ —)■ Dtt^ 
candidates, is quantihed by differences between the logarithm of likelihoods, InT/^, nnder 
hve separate mass hypotheses, h G {e, p., vr, A,p} (DLL). For the danghters from the D 
candidate, the kaon mnst satisfy DLL^^^r = ffiTx — InTj^ > 2, while the charged pion 
is reqnired to satisfy DLLj^^ < —2. Candidates with a bachelor having DLL^tt > 4 are 
selected into the B^ —)■ DK^ sample (they are said to have passed the PID requirement) 
while those that do not are placed in the B^ —)■ Dtt^ sample (they are said to have failed 
the PID requirement). 

Additional restrictions are imposed after the BDT and the PID requirements in order 
to remove specihc sonrces of backgronnd. Contribntions from gennine B^ decays that 
do not incinde a D meson are snppressed throngh a selection reqnirement on the flight 
distance signihcance, FD^, dehned as the distance between the B^ and D candidate 
vertices, divided by the uncertainty on this measurement. A reqnirement of FD^ > 2 is 
applied. The total branching fractions of B^ to four-body charmless states with a 7 r° are 
currently unmeasured and their contribution is estimated by studying the contamination 
of three-body charmless modes to the B^ —)■ spectra and scaling it 

according to the known branching fractions. The efficiency of the FD^i reqnirement is 
evalnated using simulated 6 -hadron decays to four-body charmless states with a neutral 
pion. The requirement is found to be 93% effective in the suppression of this background, 
a valne compatible with that seen in data for the three-body charmless states. From 
these stndies, it is determined that the charmless backgronnds contribnte 4 ± 1, 1 ± 1, 
4 ± 1 and 3 ± 1 candidates to the snmmed-by-charge selections of B^ —)■ [A^ 7 r^ 7 r°]oh^, 
5^—)- [7r^A=*"7r°]Dh^, 5^—)- [7r+7r“7r°]Dh^ and )■ [A+A“ 7 r°]£)/i^, respectively. 

The snppressed B^ —)■ [Tr^A^TT^joh^ decays are snbject to potential contamination 
from B^ —)■ [ 7 r+ 7 r“ 7 r°]£)/i^ and B^ —)■ [A+A^vr'^joh^ decays where one of the charged pions 
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or kaons from the D candidate is misidentified as a charged kaon or pion, respectively. 
Studies performed using simulated events demonstrate that such contamination is minimal, 
contributing 1 ± 1 candidate to each decay mode. Similarly, there 

is potential cross-feed from favoured 5^—)- [K^'K^Ti^]£,h^ decays in the suppressed ADS 
signal samples in which a and tt^ are doubly misidentihed as a tt^ and a 
respectively. This contamination is reduced by vetoing any suppressed candidate whose 
reconstructed D mass, under the exchange of mass hypotheses between the daughter kaon 
and charged pion, lies within ±30 MeV/c^ of the nominal D mass. Study of the cross-feed 
contamination in the mass sidebands of the D candidates allows for an estimate of the 
residual contamination in the signal region. After all selection requirements, this residual 
cross-feed is estimated to be (3.1 ± 0.2) x 10“^ of the total favoured )■ 
events. 

For each event, only one candidate is selected for analysis. In the 3.8% of cases where 
more than one candidate is present in an event, a choice is made by selecting the candidate 
with the decay vertex with the smallest degree of freedom. 


5 Invariant mass fit 


The observables of interest are determined with a binned maximum-likelihood fit to the 
invariant mass of the selected candidates. A total of sixteen disjoint subsamples (the 
favoured and suppressed ADS modes and the two quasi-GLW modes, separated according 
to the charge of the bachelor meson, and by the bachelor PID requirement) are htted 
simultaneously. The total probability density function (PDF) used in the ht is built from 
hve main sources, described below, representing different categories of candidates in each 
subsample. 

The B^ —)■ Dtt^ signal events are modelled through the use of a modihed Gaussian 
function. 


/(m) cx exp 


( \ 
V2cr2 ± {m- nyaL,Rj ' 


(15) 


This expression describes an asymmetric peak of mean /x and width a where the values of 
aiim < jj) and aji{m > /i) parameterise the tails of the distribution to the left and to 
the right of the peak, respectively. These signal events originate from subsamples that 
fail the bachelor PID requirement for charged kaons. Genuine B^ —)■ candidates 

that pass the PID requirement are reconstructed as )■ DK^. Since these candidates 
are reconstructed under an incorrect mass hypothesis, they represent a displaced mass 
peak with a tail that extends to higher invariant mass. Such misidentihed candidates 
are modelled by the sum of two Gaussian functions, modihed to include tail components 
similar to that of Eq. IT The two modihed Gaussian functions share a mean, but have 
two separate width parameters that are permitted to hoat. For the signal peaks, all of 
the parameters are permitted to vary, with the exception of the lower-mass tail, which is 
hxed to the value found in simulation, to ensure ht stability, and later considered as a 
source of systematic uncertainty. The same shape is used for B~ and B~^ decays, although 
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the means are allowed to be different. In addition, while the —)■ and 

)■ [ 7 r^i^^ 7 r°]£)h^ signal shapes share the same width, this parameter is permitted to 
vary for the 5^—)• [K'^K~n^]DK^ and ?• modes. 

The DK^ signal events, from the subsamples that pass the PID requirement on 

the bachelor, are modelled using the same modified Gaussian function of Eq. [T^ All of 
the shape parameters are identical to those of the —)■ modes, except for the width, 

which is fixed at (95 ± 2)% of that of the B^ —)■ Dti^ modes, based upon studies made 
using simulated events. Genuine B^ — DK^ candidates that fail the PID selection (and 
thus represent misidentified B^ — Dti^ events) are described using a fixed shape from 
simulation that is later varied to assign a systematic uncertainty. 

Partially reconstructed 6 -hadron decays are found in the invariant mass region below 
the B^ mass. However, a portion may enter the signal region. Of particular concern are 
B^ ( B °) decays involving a neutral (charged) D* meson, where the D* decays to a D ° 
and a neutral (charged) pion with this latter particle missed in reconstruction, leading 
to the same final state as in the channels of interest. The D may also decay via the 
D D *^7 channel. When the 7 is missed in reconstruction, such decays may also 

mimic the desired signal candidates. There are also further contributions from ( B °) 
decays to D ° and a neutral (charged) p or A*, where the vector meson decays into an 
(h^TT^) state from which the tt^ (tt^) is missed in reconstruction. These partially 
reconstructed decays are described by parabolic functions representative of the decays 
in question, that have been convolved with a double Gaussian to account for detector 
resolution. The yields of these background components vary independently in the ht, with 
no assumption of CP symmetry. Additionally, partially reconstructed Bg^ —)■ DK^ti^ 
decays and their charge-conjugated modes are considered as background sources to the 
ADS )■ Dh^ modes. PDFs for this background are determined from simulation and 
fixed in the invariant mass fit. The Bg^ yields are permitted to float, but CP symmetry is 
assumed given the limited interference effects due to Gabibbo suppression. 

Wrongly reconstructed D meson decays are a source of background under the signal 
peaks. These are primarily decays where the 7 r° candidate is not a daughter of the D 
meson, but is wrongly assigned as such. In the hnal hts, these contributions are modelled 
using a modified Gaussian function with a tail parameter, where this component and the 
width are permitted to vary, but the mean is hxed based on a study in data. In this study, a 
binned-maximum likelihood fit is performed to the D mass distribution in a region of ±250 
MeV/c^ about the nominal D mass [^|, where the signal and background contributions 
are modelled separately. The sPlot method 36 is used to assign signal and background 
weights to the candidates and the B^ invariant mass distribution is then plotted using the 
background weights in order to ascertain how the wrongly reconstructed D background 
contribution distributes itself in the B^ mass spectrum. This study indicates that the 
background can be described by using the function of Eq. 15 with a single tail parameter. 
As such, in the hnal hts, the wrongly reconstructed D meson background is modelled as a 
fully boating modihed Gaussian function, except for the mean that is hxed. The value of 
the hxed parameter is varied in order to assess a systematic uncertainty. 
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A linear approximation is adeqnate to describe the distribntion of combinatorial 
background across the relevant invariant mass spectrum. All —)■ DK^ modes and all 

modes share the same respective shapes, though yields vary independently. 
This allows for greater fit stability as the low statistics modes share fit information from 
the higher statistics modes. 

The measured signal yields allow the fit to determine the observables of interest. For 
example, the relationship between the yield of the decay B^ —)■ 

and the physics observables is given by 


DK+ _ 


^Dtt jdK'K'k^ jdK'K'k^ 

■ ^K/tt ' -^ADS(iC) 


^exp 


1 + 


^0 


-|_ A /CttttO 

/ ^ADS(X) 


l~l~ -^Prod . 

1 -^Prod 1 -^det 


(16) 


where eexp represents experimental selection efficiency effects and Adet are detector-related 
asymmetries (both of these are further discussed in Sect. and is the total yield of 

5^—>■ decays. In the fit, an analogous expression to Eq. 16 is used for the 

corresponding B~ decay as well as comparable equations for the other decay modes and 
their associated CP observables. 

The fit is performed such that all of the observables defined by Eqs. [Sl [To|and 13 

are free parameters. The signal yields for the decay modes of interest are presented in 
Table The uncertainties are statistical only; the systematic uncertainties are discussed 
in Sect. The corresponding invariant mass spectra, separated by the charge of the B^ 
candidate, are presented in Figs. [3 m H and m 


Table 1: The final signal yields, split in categories based on the charges of the B hadron (only 
statistical uncertainties are shown). 


B decay channel 

Yield 

B~ —)■ [A“7r+7r^]£)7r“ 

18 854 ± 176 

B- ^ 

1478± 

39 

B~ —)■ [k~ 

63 ± 

13 

B- [7r-A:+7r°]DA:- 

16 ± 

9 

B~ —>■ [7r+7r“7r°]D7r“ 

1716± 

55 

B~ —)■ [7r+7r“7r°]£)A'“ 

139 ± 

19 

B- [K+K-'k^]dti- 

509 ± 

34 

B- ^ [K+K-'K^]r)K- 

49 ± 

12 


B^ decay channel 

Yield 

B~^ —)■ [A'+7r“7r^]£)7r+ 

18 882 ± 176 

5+ ^ [A:+7r-7r0]^A:+ 

1442 ± 

39 

B^ —)■ [ti^ K~ 

25 ± 

13 

5+ ^ [ti+K-'k^]dK+ 

24 ± 

9 

B^ —)■ 

1772 ± 

55 

B^ —)■ 

125 ± 

19 

B+ [A:+A:-7r%7r+ 

541 ± 

34 

B+ [K+K-'k^]dK+ 

27 ± 

12 


6 Systematic uncertainties and results 

In addition to the sources of systematic uncertainties originating from fixed PDF parameters 
in the ht, there are several other sources that are considered. In the favoured and suppressed 
ADS modes, the ratio is hxed at 7.74% based on the measurement performed in 
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Figure 1: Invariant mass distributions of selected —)• [K^TT^TT^]£)h^ candidates, separated 

by B hadron charge. —)■ DK^ signal events are in the upper plots and B^ —)> Dtt^ events 

are in the lower plots. The solid dark (red) curve represents B^ —)• DK^ events and the solid 
light (green) curve represents B^ —)• Dn^ events. The solid (grey) shape indicates partially 
reconstructed B^ decays and the heavy dotted (red) curve indicates wrongly reconstructed D 
decays. The solid (blue) line represents the total PDF and includes the combinatorial component. 


Ref. and is assigned a systematic uncertainty of 0.22%, as per the uncertainties of that 
analysis. The DK^ versus ratio, however, is permitted to vary in the 

t [7r~^n~7r^]]:)h'^ and )■ [K~^K~7r^]]:)h'^ analyses as it must be measured for each 
mode in order to determine the RqCLW observables. 

The proportion of B^ — Dh^ samples passing or failing the PID requirements is 
determined from a sample of more than 100 million D*^ decays reconstructed as D*^ —)■ 
D7i^,D —)■ This reconstruction is performed entirely using kinematic variables 

and provides a high-purity calibration sample of and tracks. The PID efficiency 
varies as a function of track momentum, pseudorapidity and detector occupancy [^. The 
average PID efficiency of the signal is determined by reweighting the calibration spectra in 
these variables to those of the candidates in the favoured ADS sample. This average PID 
efficiency is evaluated to be 84.5% and 96.3% for kaons and pious, respectively. Systematic 
uncertainties of 0.5% and 0.8% for bachelor pions and bachelor kaons, respectively, are 
assigned to the efficiencies, which arise from the reweighting procedure. 

Due to differences in interactions with the detector material, a small negative asymmetry 
is expected in the detection of K~ and mesons. An asymmetry for pions may also 
be present and is assigned a value of (0.0 ± 0.3)% [^. The difference between the kaon 
and pion detection asymmetries is taken to be —(1.1 ± 0.4)% from studies performed in 
Ref. (^. These asymmetry values also account for the physical asymmetry of the left 
and right sides of the detector, after summing the data sets from both magnet polarities. 
There is no systematic uncertainty associated with the possible difference in number of B~ 
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Figure 2: Invariant mass distributions of selected —)• candidates. See the 

caption of Fig. for a full description. The lightly dotted (blue) line represents the combi¬ 
natorial component and the long-dashed (magenta) line indicates contributions from partially 
reconstructed —)• decays where the pion is not reconstructed. 
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Figure 3: Invariant mass distributions of selected B^ —)• [7r~^7r~7r^]Dh'^ candidates. See the 
caption of Fig. for a full description. The lightly dotted (blue) line represents the combinatorial 
component. 


and mesons, since the prodnction asymmetry Aprod is a variable parameter in the fit. 

The measured observables in the analysis are related to the ratio of relative efficiencies 
between the DK^ and B^ —modes, eB^DK/^B^D-K, independent of PID 

effects. These ratios relate the efficiency differences due to trigger, reconstruction and 
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Figure 4: Invariant mass distributions of selected —?■ [K'^K~ candidates. See the 

caption of Fig. for a full description. The lightly dotted (blue) line represents the combinatorial 
component. 


selection effects. They are measured in simulation to be (97.5 ± 3.4)% for the —)■ 

and 5^—)- [n^K^Ti^]Dh^ modes, (95.7 ±2.8)% for the 5^—)■ 
modes and (98.9±2.8)% for the 5^—)- modes. The uncertainties listed are 

based on the finite size of the simulated samples and account for the imperfect modelling 
of pion and kaon absorption rates in the detector material. 

In order to estimate the systematic uncertainties from the sources described in this 
section and in Sect. the fit is performed many times, varying each source by its assigned 
uncertainty, under the assumption that the uncertainty is Gaussian distributed. The 
spread (RMS) in the distribution of the htted value of the observables is taken as the 
systematic uncertainty. These uncertainties are summarised in Table 

The values for the coherence factor, average strong-phase differences and CP-even 
fraction reported in Refs. 15 and 18 assume a uniform acceptance across the three-body 
phase space of the D decay, which is not the case in this analysis. Studies are performed 
with amplitude models for the decays of interest and a modelling of the acceptance function 
derived from simulation to assess the impact upon these parameters arising from this 
source. It is found that in all cases the biases are negligible compared to the assigned 
uncertainties. 
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Table 2: Systematic uncertainties on the observables, multiplied by a factor of 10^. ‘PID’ refers 
to the fixed PID efficiency attributed to the bachelor tracks. ‘PDFs’ refers to the uncertainties 
based on fixed parameters in the PDF shapes that are used in the invariant mass fit. ‘Sim’ refers 
to the use of simulation to calculate relative efficiencies between the —)■ DK^ and —)■ Dtt^ 
modes, in addition to the estimated charmless background contributions and the fixed DK to Dn 
ratio on the ADS modes. ‘Anstr’ refers to the interaction and detection asymmetries. The ‘Total’ 
column represents the sum in quadrature of all of the categories of systematic uncertainties. 


A 

^ADS(K) 

A Knn^ 
"^ADS(7r) 
aKKtt^ 
^qGLW(K) 

^qGLWfA") 

AKK7r° 

^qGLW(7r) 

ATTTTTT^ 

^qGLWfvr) 
A Knn^ 

p/f-TTTr® 

^ADS(K) 

jdK'wk^ 

-^ADSfvr) 

dKKitO 

-'^qGLW 

TAnnn^ 

-'^qGLW 

^Prod 


PID PDFs Sim Ainstr Total 


3.4 

1.6 

5.1 
0.9 
0.8 
0.3 
0.4 
0.3 

0.02 

23.8 

8.1 
0.3 


39.6 
7.5 

10.2 

7.9 

2.2 

0.9 

0.9 

2.0 

0.05 

24.9 

20.7 
0.3 


8.7 

4.5 

18.8 

7.3 

1.2 
0.7 

1.4 
0.6 

0.02 

36.5 

42.5 
0.5 


5.7 
6.9 

2.1 
0.9 

4.4 

4.2 

4.2 
0.1 

0.01 

7.7 

5.3 
5.0 


41.1 

11.3 

22.1 
10.8 

5.1 

4.4 
4.6 

2.1 
0.06 
50.8 

48.3 
5.0 


The results for the observables, as determined by the £t, are 

= -0.20 ± 0.27 ± 0.04 
^ADSW = 0-438 ± 0.190 ± 0.011 
<g1ww = 0.30 ±0.20 ±0.02 
^IcLWiK) = 0.054 ± 0.091 ± 0.011 
= -0.030 ± 0.040 ± 0.005 
= -0.016 ± 0.020 ± 0.004 

^ Q Q^Q ^ Q Q2g ^ Q QQ5 

RadsIk) = 0.0140 ± 0.0047 ± 0.0021 
-Radsw = 0.00235 ± 0.00049 ± 0.00006 
= 0.95 ± 0.22 ± 0.05 
= 0.98 ±0.11 ±0.05 
-^Prod -0.0008 ± 0.0055 ± 0.0050 
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where the hrst uncertainties are statistical and the second are systematic. 

None of the asymmetry observables exhibit any signihcant CP violation. The results for 


1 TTTTTT 

^qGLW{if) 


this is the hrst time 


the ADS observables are more precise than those obtained by previous experiments 19,20 
and are compatible with them. Furthermore, apart from 
that the quasi-GLW observables have been measured. 

A likelihood-ratio test is used to assess the signihcance of the suppressed ADS signal 
yields, as well as those of the —)■ [K^K~ decays. This is performed by calculating 

the quantity ■>/—2 ln(£b/-^^s+b) where £b and £s+b are the maximum likelihood values of the 
background-only and signal-plus-background hypotheses, respectively. Including systematic 
uncertainties, signihcances of 5.3a and 2.8a are found for the —>■ [ 7 r^iF^ 7 r°]£) 7 r^ and 

—)■ decays, respectively. For the —)■ [K^K~7T^]uh'^ selections, the 

—>■ mode is found to have a signihcance greater than lOa, while a signihcance of 

4.5a is measured for the 5^—> decay. 


7 Interpretation and conclusions 



20 40 60 80 100 120 140 160 180 

Y [degrees] 

Figure 5: Scan of the probabilities over the parameter space. Shown are the na profile 
likelihood contours, where = n?, with n = 1 being the light (blue) shaded region, n = 2 the 
dark (blue) shaded region and n = 3 corresponding to the white area. The result is seen to be 
compatible with the current LHCb measurement of 7 and rs, indicated by the point with error 
bars. 

The measured observables from the B^ —)■ DK^ decay channels are used to obtain 
constraints on the underlying physics parameters r^, Sb and 7 . For this purpose, the small 
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Figure 6 : Scan of the probabilities over the 7-^5 parameter space. Shown are the na profile 
likelihood contours, where Ax^ = n?, with n = 1 being the light (blue) shaded region, n = 2 the 
dark (blue) shaded region and n = 3 corresponding to the white area. The result is compatible 
with the current LHCb measurement of 7 and 5b, indicated by the point with error bars. 


effects of mixing and interference in —)■ decays are neglected. Using the 

measurements and associated fit covariance matrix and systematic uncerta inty correlations, 

?K+K-tt° 


and taking external measurements of 


15 


18 


and the 


and 

branching ratios of the D decay channels as additional inputs with their associated 
uncertainties, a global minimisation is performed. A scan of the physics parameters 
is executed for a range of values and the difference in goodness of ht, between the 
parameter scan values and the global minimum, is evaluated. Assuming that this 
minimisation function is distributed in a Gaussian manner enables a probability to be 
assigned for each set of values of the physics parameters. 

Two-dimensional scans are performed for 7 vs. tb and 7 vs. 6b in the ranges 0.03 < 
tb < 0.16, 0° < 5b < 180° and 0° < 7 < 180°. Figs. and shows the la, 2a and 3a 
contours determined from these scans. It can be seen that the results are compatible with 
the values obtained from a global analysis of other LHCb measurements sensitive to 7 at 
tree level , which are also shown. The scans return a best-ht value for the parameter 
Tb of 0.11 ± 0.03. No useful constraints are obtained for either 7 or 6b- However, the 
measurements of the observables are expected to provide improved precision on these 
parameters when included in a global analysis of all LHCb —)■ DK^ results. 

In summary, measurements of CP asymmetries and related observables have been 
performed using B^ —)■ DK^ and B^ —)■ decays with an inclusive analysis of the ADS 

modes D —and, for the hrst time, the quasi-GLW modes D —)■ 7 r+ 7 r“ 7 r° and 
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D —)■ The results for the ADS observables are the most precise measurements 

of these quantities. No evidence of CP violation is obtained with the current experimental 
precision. First observations have been made of the decays —)■ [ 7 r^iF^ 7 r°]£) 7 r^ and 

[K^K~, and hrst evidence is obtained for the mode —?• [K^. 

When analysed in the context of the underlying physics parameters, the results exhibit 
good consistency with other LHCb measurements. The measurements will be valuable in 
improving knowledge of the unitarity triangle angle 7 when combined with LHCb results 
from —)■ DK^ measurements exploiting other D decay channels. 
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